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O . Abstract 

\o ■ 

^^ I In this paper we consider the matching coefficients up to two loops between 

Quantum Chromodynamics (QCD) and Non-Relativistic QCD (NRQCD) for the 
vector, axial- vector, scalar and pseudo-scalar currents. The structure of the effective 



Oh, theory is discussed and analytical results are presented. Particular emphasis is put 

on the singlet diagrams. 



PACS numbers: 

In the recent years quite a lot of activity has been devoted to the treatment of bound 
states of two heavy particles both in QED and QCD (for a recent review see, e.g., Ref. [1]). 
From the theory point of view the calculations have been put onto a solid basis due to the 
formulation of proper effective theories [2,3], NRQED and NRQCD, respectively, which 
provide the possibility to systematically evaluate higher order corrections. The construc- 
tion of the effective theories consists of essentially two steps: First, the effective operators 
involving the light degrees of freedom have to be constructed and second the correspond- 
ing couplings, the so-called coefficient functions, have to be computed by comparing the 
full and the effective theories. The latter is also referred to as matching calculation. 

The framework which is considered in this letter consists of QCD accompanied by ex- 
ternal currents where we allow for vector, axial- vector, scalar and pseudo-scalar couplings. 
The main results of this letter are the two-loop matching coefficients. Thus, following the 
prescription outlined above we determine in a first step the effective currents and then 
perform a matching calculation. 

The matching coefficients provided in this paper constitute a building block in all 
calculations involving the corresponding external currents. This includes in particular 
production and decay processes of heavy quarkonia or the production of top quark pairs 



close to threshold. One could also think of the decay of a CP-even or CP-odd Higgs boson 
(with mass M) into two quarks with 2m ^ M. 

The basic idea behind the construction of the Lagrange density for NRQCD is to 
expand all terms of the QCD Lagrangian in the limit of a large quark mass. A similar 
procedure has to be applied to external currents which we define in coordinate space as 
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= ^7^^, 
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= V^7^75^ 


3s 


= V^V", 


Jp 


= ^«75^ • 



(1) 

Note that the anomalous dimension of j^ and j^ is zero whereas for the scalar and pseudo- 
scalar current it is obtained from the renormalization constant Zg = Zp = Z^, which is 
given at the two- loop level in Ref. [4]. 

In order to perform the transition to the effective theory it is convenient to work in 
momentum space and to introduce the two-component Pauli-spinors in the form 



where m denotes the heavy quark mass. In Eq. Q x is a spinor that annihilates a heavy 
quark and correspondingly creates a heavy anti-quark with momentum p. 

In a first step we want to express the currents of Eq. (Q) in terms of and x ^ind 
expand in the inverse heavy quark mass. This actually leads to the tree-level matching 
conditions. When inserting Eq. (jSJ into (P) it turns out to be convenient to split the 
time-like and space-like coefficients of the vector and axial- vector currents. This leads to 
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where k = 1,2,3 and the currents in the effective theory are given by 
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= ^J^[a^p-a]x, 
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(4) 

Note that jp also appears in the expansion of j^ which means that the corresponding 
matching coefficients are equal. This will be used as a check of our calculation. Due to 
the occurrence of the momentum p in j^ and js an expansion in the external momenta 
has to be performed in order to obtain the loop corrections to the corresponding matching 
coefficients. 

The basic idea to obtain the matching coefficients is to compute vertex corrections 
induced by the considered current both in the full and the effective theory. In practice 
it is convenient to consider the renormalized vertex function with two external on-shell 
quarks and to perform an asymptotic expansion about s = 4m^, where s is the momentum 
squared of the external current, the so-called threshold expansion [5,6]. Denoting by T^ the 
proper structure of the genuine vertex corrections and by Z2 and Z^ the renormalization 
constants due to the quark wave function and the anomalous dimension of the current 
one obtains the equation 

Z2Z^T^{qi, ga) = c^Z2Z;^^f^ + ... , (5) 

where x G {v, a, s,p} with the understanding that the axial-vector part is split into time- 
like and space-like components. The ellipses denote terms suppressed by inverse powers 
of the heavy quark mass and the quantities in the effective theory are marked by a tilde. 
Cx is the matching coefficient we are after. In our approximation Z2 = 1. Z2 to two 
loops has been computed in Ref. [7]. As far as T^ is concerned only the tree- level result 
determined by jx contributes to Eq. (0). The momenta gi and g2 in Eq- ® correspond 
to the outgoing momenta of the quark and anti-quark which are considered on-shell. 

Starting from order a^ the matching coefficients c^ exhibit infra-red divergences which 
are compensated by ultra-violet divergences of the effective theory rendering physical 
quantities finite. In Eq. (jSJ the renormalization constant Z^ which generates the anoma- 
lous dimension oi jx takes over this part. 

The quantities Tx are conveniently obtained with the help of projectors which are 
constructed in such a way that they project on the coefficients of Tx. For the vector 
case, the zeroth component of the axial-vector and the pseudo-scalar case we can simply 
identify qf = q2 = ?^/4 = m^ and use 

r„ = Tr [p^^)r(^)''^] , 

Tp = Tr [p(p)r(P)] , 
r,,o = Tr [P^-'0)r(")'^] , (6) 



Figure 1: Feynman diagrams contributing to the matching coefficients. In (e) the so- 
called singlet diagram is shown which does not contribute to c^. In the closed fermion 
loop all quark flavours have to be considered. 
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As already mentioned above the case (a, 0) is used as a check for the pseudo-scalar match- 
ing coefficient. 

For the axial- vector and scalar cases we have the equations analogous to Eq. ©. 
However, since the corresponding effective currents have a suppression factor |p |/m it is 
necessary to choose gi = g/2 + p and g2 = Q'/2 — p, to expand up to linear order in p 
and to set afterwards p = and q^ = Am?. Note that we choose a reference frame where 
g -p = [5,6]. Thus the projectors are more complicated and are given by 
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In Fig. ^ some Feynman diagrams contributing to the matching coefficients are shown. 
Due to the application of the projectors the corresponding integrals can be reduced to 
the functions J± and L± as defined in Eqs. (14) and (55) of Ref. [5]. However, due 
to the expansion in the momentum p the powers of the denominators are higher and 
a systematic reduction of the scalar integrals to master integrals is necessary. For the 
current calculation we implemented the method of Ref. [8]. For some of the occuring 
integrals the program AIR [9] is applied. The details will be described elsewhere. 

An important class of diagrams is constituted by the so-called singlet diagrams (cf. 
Fig. n^e)) where the external current does not couple to the quark-anti-quark pair of 



the final state. Due to Furry's theorem there is no contribution to the vector case from 
these diagrams, however, non-vanishing, finite results are obtained for c^, c^ and Cp. For 
the scalar and the pseudo-scalar currents only the heavy quark is running in the closed 
fermion loop. All other quark flavours are suppressed by the light quark mass. This is 
different for the axial-vector coupling. Here we consider the effective current formed by 
the top and bottom quark field 

Ja^ = i^^l^t-ln^h, (9) 

which ensures the cancellation of the anomaly-like contributions. For the same reason the 
contributions from the remaining light quarks cancel. 

In the analytical results given below the contributions from the singlet diagrams are 
marked separately. At this point we only want to mention that in the axial-vector and 
pseudo-scalar case 75 was treated according to the prescription of Ref. [10]. In practice 
this means that we perform the replacements 

75 - ^e'^^^-7M7.7p7. , (10) 

strip off the e tensor and deal with the objects with three and four indices, respectively. 
The corresponding projectors are obtained by performing the replacements of Eq. (jlOj) in 
Eqs. (|7j) and (jH)) which makes them more complicated. However, the very calculation is 
in close analogy to the non-singlet case. After summing all two-loop contributions one 
obtains a finite result. 

Note that for the non-singlet contributions it is save to use anti-commuting 75. Ac- 
tually, the treatment according to Ref. [10] leads to a wrong result. This is due to the 
infra-red divergences which are absent in the singlet diagrams. 

An alternative method to perform the calculation of the vertex corrections is based on 
the evaluation of the tensor integrals. Here, we used the T-operator method of Ref. [11] 
to reduce the tensor integrals to products of the metric tensor and external momenta and 
scalar integrals with shifted space-time dimension. The resulting Dirac structures were 
further simplified and rewritten in terms of NRQCD fields. 

In addition to the bare two-loop diagrams we have to take into account the one-loop 
renormalization contribution from the heavy quark mass, which we renormalize on-shell, 
and the strong coupling renormahzed in the MS scheme. 

We want to mention that all contributions have been evaluated for general gauge 
parameter ^. The final results for the matching coefficients are all independent of ^ which 
constitutes an important check on our calculation. 

Let us in the following present our results and compare with the literature. The 
two-loop matching coefficient for the vector current has been computed almost ten years 
ago [12, 13]. We confirmed these results and provide for completeness the analytical 



expressions 
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where Ca = Nc and Cp = {N^ — l)/{2Nc) are the Casimir operators of the adjoint and 
fundamental representation of SU(A^c), respectively, T = 1/2, and rii is the number of 
massless quarks, (n denotes Riemann's zeta-function. The one-loop result can already 
be found in Ref. [14]. The anomalous dimension of the effective vector current, which is 
related to Z^ through 7^, = ^^j^, reads 



7v = -(^)^(2C^ + 3C^Ca)C2. 

Our results for the two-loop matching coefficients c^, Cg and Cp are given by 
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The one-loop result for Cp can already be found in Ref. [15]; the two-loop coefficients of 
Eq. (|T3| are new. They constitute our main result. The one- loop coefficients can be easily 
obtained from the one-loop on-shell vertex corrections with arbitrary momentum squared 
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of the external current, s. In the analytic expressions it is straightforward to perform the 
limit where the velocity of the produced quarks is small. After subtracting the leading 
term, which corresponds to the Coulomb singularity, one remains with the result for the 
matching coefficients [16]. At two loops this simple trick does not work any more and the 
calculation has to be performed from scratch as has been done in this letter. 
The contributions from the singlet diagrams correspond to 

^s^, = -^C2 + 4C2ln2-21n2 + ^ln22 + m('l-^ln2 

41 = ---C2 + 4C2ln2-ln2 + z-, 

^fng = ^C2 + 3C2ln2-|c3 + ^^^C2. (14) 

Xgi^ and Xgi^ receive only contributions from diagrams which are finite and contain only 
the heavy quark. The corresponding result holds both for top and bottom quarks. This 
is different in the case of X^j^ . Actually, the result in Eq. (fT^ corresponds to the case 



where top quarks are considered in the final state. Note that X^^' receives contributions 



(a) 



sing 



from diagrams with top and bottom quarks in the closed triangle loop (cf. Fig. life)). 
Taken separately they are divergent, however, the sum is finite. If one considers bottom 
quarks in the final state one still has to consider top and bottom quarks in the closed 
triangle loop. Again only the sum of all diagrams is finite with the result 

The diagram in Fig. ^e) for axial- vector coupling and with bottom in the final state was 
also considered in Ref. [17, 18], for abritrary values of s and m^, but for rrib = so that 
no direct comparison with Eq. (jl5|) is possible. 

As mentioned above, it is possible to extract the result for Cp from the zero component 
of the axial- vector current. This is quite evident in the non-singlet case. However, for the 
singlet contribution this check is highly non-trivial since in this approach Cp is obtained 
from diagrams both with top and bottom quarks in the closed triangle whereas in the 
direct calculation only one type of quarks appears. 

The singlet results of Eqs. fll4|l and ()15|) and the fermionic contributions of Eq. p3j) 
are in agreement with Ref. [19-21] where the off-shell contributions have been considered. 
Since they do not develop an infrared singularity the limit s -^ Am? can be performed. 
This is different in the case of the non-fermionic contributions where due to the infrared 
divergence the off-shell results [20-22] cannot be used in order to obtain the matching 
coefficients. 

For completeness we also provide the result for the anomalous dimensions correspond- 



ing to Eq. (jT^ which read 

Is = -(^)'(4C^ + C^C^)C2, 

7p = -(^)'(6C^ + 3C^C^)C2. (16) 

The result for 7p agrees with the one extracted from Ref. [23]. 

We want to mention that the coefficient Cp has been considered in Ref. [24] in the 
context of the Be meson. The latter consists of two heavy quarks, however, with different 
masses. This makes the calculation significantly more difficult since two instead of one 
mass scale appear in the integrals. In Ref. [24] the reduction to master integrals has been 
performed exactly whereas the latter have been evaluated in the limit rric <^ rrih so that 
a comparison with the present analysis is not possible. 

To summarize, in this paper we computed the two-loop matching coefficients between 
QCD and NRQCD for an axial-vector, scalar and pseudo-scalar current. Furthermore, we 
performed an independent check of the matching coefficient in the vector case. The latter 
contributes to the second order result of the threshold production of top quark pairs. The 
result for the axial-vector current only contributes to the fourth-order analysis which is 
currently still out of reach. 
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